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Abstract 



We study integrable lattice regularizations of the Sine-Gordon model with the help of 
the Separation of Variables method of Sklyanin and the Baxter Q-operators. This leads 
us to the complete characterization of the spectrum (eigenvalues and eigenstates), in 
terms of the solutions to the Bethe ansatz equations. The completeness of the set 
of states that can be constructed from the solutions to the Bethe ansatz equations is 
proven by our approach. 
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1. Introduction 
1.1 Motivation 

The study of the Sine-Gordon model has a long history. It has in particular served as an im- 
portant toy model for interacting quantum field theories. The integrability of this model gives 
access to detailed non-perturbative information about various characteristic quantities, which 
allows one to check physical ideas about quantum field theory against exact quantitative results. 

It is particularly fascinating to compare the Sine-Gordon model with the Sinh-Gordon model. 
The Hamiltonian density hsa of the Sine-Gordon model and the corresponding object hshc of 
the Sinh-Gordon model, 



are related by analytic continuation w.r.t. the parameter /3 and setting /3 = ib. The integrability 
of both models is governed by the same algebraic structure Uq{sl2) with q = e~^^^^ . This leads 
one to expect that both models should be closely related, or at least have the same "degree of 
complexity". 

The physics of these two models turns out to be very different, though. Many of the key objects 
characteristic for the respective quantum field theories are not related by analytic continuation 
in the usual sense. While the Sine-Gordon model has much richer spectrum of excitations and 
scattering theory in the infrared (infinite R) limit, one may observe rather intricate structures 
in the UV- limit of the Sinh-Gordon model [|Za06ll . which turn out to be related to the Liouville 
theory [|ZZ95[ ITOSai IBT091 . These differences can be traced back to the fact that the periodicity 
of the interaction term Svr/^ cos(2/90) of the Sine-Gordon model allows one to treat the variable 
as angular variable parameterizing a compact space, while is truly non-compact in the 
Sinh-Gordon model. 

The qualitative differences between the Sine-Gordon and the Sinh-Gordon model can be seen 
as a simple model for the differences between Nonlinear Sigma-Models on compact and non- 
compact spaces respectively. This forms part of our motivation to revisit the Sine-Gordon model 
in a way that makes comparison with the Sinh-Gordon model easier. 

1.2 Open problems 




h 



h 



ShG — 



+ {d^(pY + Svr/i cos(2/30) , 
+ {d^(pf + 87r/i cosh(260) , 



(1.1) 



A lot of important exact results are known about the Sine-Gordon model. Well-understood are 
in particular the scattering theory in the infinite volume. The spectrum of elementary particle 
excitations and the S-matrix of the theory are known exactly [|KT77[ IZa77[ IFSTSOi IK08OL 
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Relatedly, there is a wealth of information on the form-factors of local fields, see e.g. [|Sm92[ 
IBFKZllLZOlU for the state of the art and further references. In the case of finite spacial volume, 
the nonlinear integral equation^] derived by Destri and De Vega [IDDV921 IDDV94[ IDDV97[ 
FMQR971 |FRT98[ rFRT99ll give a powerful tool for the study of the finite-size corrections to the 
spectrum of the Sine-Gordon model. 

However, there are several questions, some of them fairly basic, where our understanding does 
not seem to be fully satisfactory. We do not have exact results on correlation functions on the 
one hand, or on expectation values of local fields in the finite volume on the other hand at 
present. 

Even the present level of understanding of the spectrum of the model does not seem to be fully 
satisfactory. The truth of the commonly accepted hypothesis that the equations derived by Destri 
and De Vega describe all of the states of the Sine-Gordon model has not been demonstrated yet. 
The approach of Destri and De Vega is based on the Bethe ansatz in the fermionized version 
of the Sine-Gordon model, the massive Thirring model nDDV87ll . This approach a priori only 
allows one to describe the states with even topological charge, and it inherits from its roots in 
the algebraic Bethe ansatz some basic difficulties like the issue of its completeness. 

In the Bethe ansatz approach it is a long-standing problem to prove that the set of states that is 
obtained in this way is complete. Early attempts to show completeness used the so-called string 
hypothesis which is hard to justify, and sometimes even incorrect. At the moment there are 
only a few examples of integrable models where the completeness of the Bethe ansatz has been 
proven, including the XXX Heisenberg model, see HMTVH and references therein. A similar 
result has not been available for the Sine-Gordon model or its lattice discretizations yet. One of 
the main results in this paper is the completeness result for the lattice Sine-Gordon model. We 
prove a one-to-one correspondence between eigenstates of the transfer matrix and the solutions 
to a system of algebraic equations of the Bethe ansatz type. For brevity, we will refer to this 
result as completeness of the Bethe ansatz. We furthermore show that the spectrum of the 
transfer matrix is simple in the case of odd number of lattice sites, and find the operator which 
resolves the possible double degeneracy of the spectrum of the transfer matrix in the case of 
even number of lattice sites. 

1.3 Our approach 

We will use a lattice regularization of the Sine-Gordon model that is different from the one used 
by Destri and De Vega. It goes back to HFSTSOi ITK82L and it has more recently been studied in 
nF94[|FV94ll . For even number of lattice sites the model is related to the Fateev-Zamolodchikov 

'This type of equations were before introduced in a different framework in ||KP91||KBP91| 
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model IIFZ82L as was observed in IIFV94II . or more generally to the Chiral Potts model, as 
discussed in the more recent works [BBR961 IBaOSH . This allows one to use some powerful 
algebraic tools developed for the study of the chiral Potts model nBS90l in the analysis of the 
lattice Sine- Gordon model. 

The issue of completeness of the Bethe ansatz had not been solved in any of these models yet. 
What allows us to address this issue is the combination of Separation of Variables method (SOV- 
method) of Sklyanin HSkSSl |Sk92[ |Sk95l with the use of the Q-operators introduced by Baxter 
[|Ba72L We will throughout be working with a certain number of inhomogeneity parameters. It 
turns out that the SOV-method works in the case of generic inhomogeneity parameters where the 
algebraic Bethe ansatz method fails. It replaces the algebraic Bethe ansatz as a tool to construct 
the eigenstates of the transfer matrix which correspond to the solutions of Bethe's equations. In 
a future publication we will show that the results of our approach are consistent with the results 
of Destri and De Vega. 

Another advantage of the lattice discretization used in this paper which may become useful in 
the future is due to the fact that one directly works with the discretized Sine-Gordon degrees of 
freedom, which is not the case in the lattice formulation used by Destri and De Vega. Working 
more directly with the Sine-Gordon degrees of freedom should in particular be useful for the 
problem to calculate expectation values of local fields. This in particular requires the deter- 
mination of the SOV-representation of local fields analogously to what has been done in the 
framework of the algebraic Bethe ansatz in liKMT99l IMTOOII . The SOV-method in principle 
offers a rather direct way to the construction of the expectation values, as illustrated in the case 
of the Sinh-Gordon model by the work ULuOlll . 

Acknowledgements. We would like to thank V. Bazhanov and F. Smirnov for stimulating discussions, 
and J.-M. Maillet for interest in our work. 

We gratefully acknowledge support from the EC by the Marie Curie Excellence Grant MEXT-CT-2006- 
042695. 

2. Definition of the model 

2.1 Classical Sine-Gordon model 

The classical counterpart of the Sine-Gordon model is a dynamical system whose degrees of 
freedom are described by the field 0(2;, t) defined for (x, t) E [0, i?] x M with periodic bound- 
ary conditions (f){x + R,t) = The dynamics of this model may be described in the 
Hamiltonian form in terms of variables (f)(x, t), n(x, t), the Poisson brackets being 

{ n(x, t) , t)} = 2n 6{x - x') . 



The time-evolution of an arbitrary observable 0(t) is then given as 
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dtO{t) = {H,0{t)}, 
with Hamiltonian H being defined in (|l.ll) . 

The equation of motion for the Sine-Gordon model can be represented as a zero curvature 
condition, 

[dt-V{x,t;X),d,-U{x,t;X)] = 0, (2.1) 
with matrices U{x, t; A) and V{x, t; A) being given by 



U{x,t;\) 
Vix,t;X) 



if n -zm(Ae-^'^'^ - A-^e*^' 

if0' +im{Xe-'^'f' + X-^e'^'^) 



+im{Xe'^''' + A-^e"*'^'^) -if 0' 

and m related to fihy = tt/S'^ij,. 



(2.2) 



2.2 Discretization and canonical quantization 

In order to regularize the ultraviolet divergences that arise in the quantization of these models 
we will pass to integrable lattice discretizations. First discretize the field variables according to 
the standard recipe 

(f)n = <P{nA) , n„ = An(r2A) , 

where A = _R/N is the lattice spacing. In the canonical quantization one would replace (pn, n„ 
by corresponding quantum operators with commutation relations 

[0„ , n„] = 2'n:i6n,m ■ (2.3) 

Planck's constant can be identified with by means of a rescaling of the fields. 

The scheme of quantization of the Sine-Gordon model considered in this paper will deviate 
from the canonical quantization by using u„ = e^2^" and v„ = e'^'^^" as basic variables. For 
technical reasons we will consider representations where both u„ and v„ have discrete spectrum. 
Let us therefore take a moment to explain why one may nevertheless expect that the resulting 
quantum theory will describe the quantum Sine-Gordon model in the continuum limit. 

First note (following the discussion in IIZa94ll ) that the periodicity of the potential 8nfi cos(2/30) 
in (11.11) implies that shifting the zero mode 0o = dx (f){x) by the amount 7r//3 is a symme- 

i Dp. 

try. In canonical quantization one could build the unitary operator \N = e^i^ which generates 
this symmetry out of the zero mode po = J^^ dx Il{x) of the conjugate momentum 11. W 
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should commute with the Hamihonian H. One may therefore diagonaUze W and H simuhane- 
ously, leading to a representation for the space of states in the form 

■H - [ da where \N ■ = e^^H^ . (2.4) 

Jsi 

An alternative way to take this symmetry into account in the construction of the quantum theory 
is to construct the quantum theory separately for each a-sector. This implies that the field 
should be treated as periodic with periodicity vr//3, and that the conjugate variables n„ have 
eigenvalues quantized in units of 13, with spectrum contained in { 2a/3/N + 47r/3fc ; k E Z,}. The 
spectrum of n„ is such that the operator W = e"^^^", with Rpo approximated by J2n=i is 
realized as the operator of multiplication by e*°. 

Let us furthermore note that it is possible, and technically useful to assume that the lattice 
field observable 0„ has discrete spectrum, which we will take to be quantized in units of 13. In 
order to see this, note that the field is not a well-defined observable due to short-distance 
singularities, whereas smeared fields like Jjdx(f){x), I C [0,-R] may be well-defined. The 
observable Jj dx would in the lattice discretization be approximated by 

~ 5^ A0„. (2.5) 

nAel 

So even if 0„ is discretized in units of (3, say, we find that the observable is quantized in 
units of A/3, which fills out a continuum for A — t- 0. 



2.3 Non-canonical quantization 

As motivated above, we will use a quantization scheme based on the quantum counterparts of 
the variables u„, n = 1, . . . , N related to n„, 0„ as 

= e'^^" , Vn = e"^^*" . (2.6) 

The quantization of the variables m„, f „ produces operators u„, which satisfy the relations 

UnVm = g''"™v^u„ , wheie q = e""^' . (2.7) 

We are looking for representations for the commutation relations (12.71) which have discrete 
spectrum both for u„ and v„. Such representations exist provided that the parameter g is a root 
of unity, 

f3' = ^, p,pE'L>\ (2.8) 
P 

We will restrict our attention to the case p odd and p' even so that = 1. It will often be 
convenient to parameterize p as 

p = 21 + 1, / e . (2.9) 
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Let us consider the subset §p = {g^"; n = 0, . . . , 21} of the unit circle. Note that E>p = {g"; n = 
0, . . . , 2/} since g^'"''^ = q. This allows us to represent the operators u^, v„ on the space of 
complex- valued functions ^ : — )■ C as 

Un ■ 1p{zi, . . . , Zn) = UnZn1p{Zi, . . . , . . . , Zn) , jq-j 
V„ ■ . . . , 2n) = Vni^iZi, q~^Zn, • • • , ^n) • 

The representation is such that the operator u„ is represented as a multiplication operator. The 
parameters Un, Vn introduced in (|2.10l) can be interpreted as "classical expectation values" of 
the operators u„ and v„. The discussion in the previous subsection suggests that the Vn will be 
irrelevant in the continuum limit, while the average value of m„ will be related to the eigenvalue 

of W via m„ = exp(z/32a/N). 

2.4 Lattice dynamics 

There is a beautiful discrete time evolution that can be defined in terms of the variables in- 
troduced above which reproduces the Sine-Gordon equation in the classical continuum limit 
nFV94|| . It is simplest in the case where m„ = 1, f„ = 1, n = 1, . . . , N. We will mostl>§ restrict 
to this case in the rest of this paper. 

More general cases were treated in [|BBR96[lBa08ll . 
2.4.1 Parameterization of the initial values 

As a convenient set of variables let us introduce the observables defined as 

— ^-2if)4>n ^ — gi^(n„+n„_i-2(/)„— 2(;i„_i) ^2 \ \^ 



These observables turn out to represent the initial data for time evolution in a particularly con- 
venient way. The quantum operators f„ which correspond to the classical observables /„ satisfy 
the algebraic relations 

^2n±lhn = 'fhnhn±l, q = e , ^n^n+m = ^n+m^n for m>2. (2.12) 

There exist simple representations of the algebra (I2.12|) which may be constructed out of the 
operators u„, v„, given by 

hn = V^, f2„_l = U„U„_i. (2.13) 

The change of variables defined in (12.131 ) is invertible if N is odd. 



^Except for Section H) 
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2.4.2 Discrete evolution law 

Let us now describe the discrete time evolution proposed by Faddeev and Volkov [IFV94II . 
Space-time is replaced by the cylindric lattice 

£ = {{u,t),u e Z/NZ , r G Z , z/ + r = even } . 

The condition that z/ + r is even means that the lattice is rhombic: The lattice points closest to 
(z/, r)are(z/±l,r + l) and (z^ ± 1, r — 1). We identify the variables f„ with the initial values 
of a discrete "field" ^y^^ as 

One may then extend the definition recursively to all (z/, r) G £ by means of the evolution law 

^u,T+i = g^{q^u~i,r) ■ K,T-i ■ 9^{q^y+i,r) , (2-14) 

with function g defined as 

= ^ (2.15) 

where k plays the role of a scale-parameter of the theory. We refer to nFV94|| for a nice dis- 
cussion of the relation between the lattice evolution equation (12.141) and the classical Hirota 
equation, explaining in particular how to recover the Sine-Gordon equation in the classical con- 
tinuum limit. 

2.4.3 Construction of the evolution operator 

In order to construct the unitary operators U that generate the time evolution (12.141) let us intro- 
duce the function 

^ -I < \ 2r—l 
r=l 

which is cyclic, i.e. defined on Zp. The function Wx{z) is a solution to the functional equation 

{z + X)Wx{qz) = (1 + Xz)Wx{q-'z) , (2.17) 
which satisfies the unitarity relation 

{Wx{z)r = iWx,{z))-\ (2.18) 

Note in particular that Wx{z) is "even", i.e. Wx{z) = Wx{l/z). Further properties of this 
function are collected in Appendix A. 



11 

Let us then consider the operator U, defined as 

N N 

U = n W^«-(f2n) ■ Uo ■ n W^«-(f2n-l) , (2.19) 
n=l n=l 

where Uq is the parity operator that acts as Ug ■ = f^T^ ■ Ug. It easily follows from (12.171) that 
U is indeed the generator of the time-evolution (12.141) . 

t,r+l = ■ ■ U . (2.20) 

One of our tasks is to exhibit the integrability of this discrete time evolution. 



3. Integrability 

The integrability of the lattice Sine-Gordon model is known [|lK82l IFV92[ IBKP93[ IBBR96L 
The most convenient way to formulate it uses the Baxter Q-operators nBa72i These operators 
have been constructed for the closely related Chiral Potts model in PBS90II . By means of the 
relation between the lattice Sine Gordon model and the Fateev-Zamolodchikov model summa- 
rized in Appendix ID] one may adapt these constructions to the formulation used in this paper. 
For the reader's convenience we will give a self-contained summary of the construction of the 
T- and Q-operators and of their relevant properties in the following section. 



3.1 T-operators 

As usual in the quantum inverse scattering method, we will represent the family Q by means 
of a Laurent-polynomial T(A) which depends on the spectral parameter A. The definition of 
operators T(A) for the models in question is standard. It is of the general form 

T(A) = trc2M(A) , M(A) = ^^(A/^n) • • • ^i(A/6) , (3.1) 

where we have introduced inhomogeneity parameters . . . , as a useful technical device. 
The Lax-matrix may be chosen as 



,-lw-l 



rSG/\^ _ I '''^nil ^/^n^n + Q'^^'^n^^n^) "^n^n "n 1 on 

[■^) — — \ . _i ,1 . _!/ ,1 _i , _l _l^ • (-'•^) 

\ A V —A V tu (7 2 K V +0 2/£ V 

n n n n n n n ' ^ n n j 



An important motivation for the definitions (|3.1I) . (|3.2I) comes from the fact that the Lax-matrix 
L^{\) reproduces the Lax-connection U (x) in the continuum limit. 

The elements of the matrix M(A) will be denoted by 

^ /a(a) b(aA 

U(A) D{\) 



They satisfy commutation relations that may be summarized in the form 
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R{X/li) (M(A) ® 1) (1 ® M(^)) = (1 ® M(/i)) (M(A) ® l)R{\/^l) , 



(3.4) 



where the auxiliary R-matrix is given by 



/ q\-q-^\ 



\ 



R{X) 



q-q 



(3.5) 



q-q 



1 A-A-i 



V 



qX~q 



It will be useful for us to regard the definition (13.11) as the construction of operators which 
generate a representation 7?.n of the so-called Yang-Baxter algebra defined by the quadratic 
relations (13.41) . The representation TZ^ is characterized by the 4N parameters n = (ki, . . . , kn), 
^ = (^1, • • • , ^n), u = {ui, . . . , mn) and = (vi, wn). 

The fact that the elements of M(A) satisfy the commutation relations (13.41) forms the basis for 
the application of the quantum inverse scattering method. The mutual commutativity of the 
T-operators, 



follows from (13.41) by standard arguments. The expansion of T(A) into powers of A produces 
N algebraically independent operators Ti, . . . , Tn. Our main objective in the following will be 
the study of the spectral problem for T(A). The importance of this spectral problem follows 
from the fact that the time-evolution operator U of the lattice Sine-Gordon model will be shown 
to commute with T(A) in the next section. 

3.2 Q-operators 

Let us now introduce the Baxter Q-operators Q(/i). These operators are mutually commuting 
for arbitrary values of the spectral parameters A and /x, and satisfy a functional relation of the 
form 



with a(A) and d{\) being certain model-dependent coefficient functions. The generator of lat- 
tice time evolution will be constructed from the specialization of the Q-operators to certain 
values of the spectral parameter A, making the integrability of the evolution manifest. 



[T(A),T(/x)] = 0, 



(3.6) 



T(A)Q(A) = a(A)Q(g-U) + d(A)Q(gA), 



(3.7) 
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3.2.1 Construction 



In order to construct the Q-operators let us introduce the following renormalized version of the 
function Wx{z), 

'"^to'") = n^;^n^B^. (3.8) 

r=l ^ r=l ^ 

The function wx{z) is the unique solution to the functional equation (12.171) which is a polyno- 
mial of order / in A and which satisfies the normalization condition wi(g^") = 1. 

The Q-operators can then be constructed in the form 

Q(A,/i) = Y(A)-(Y(/i*))t, (3.9) 
where Y(A) is defined by its matrix elements Y^(z, z') = ( z | Y(A) | z' ) which read 

N 
n=l 

where e = —iq~^, and wx{z) is the discrete Fourier transformation of w{z), 

^ I I 
^^(^) = - XI ^''^^(^'"^ ' = ^y'^^wxiq"^). (3.11) 

^ r=-l r=-l 

Note in particular the normalization condition Wi{q'^) = Sr^. 

Despite the fact that Q(A,/i) is symmetric in A and fi, Q(\, fi) = Q(yU, A) as follows from the 
identity (IB. 61) proven in Appendix B, we will mostly consider as a fixed parameter which will 
later be chosen conveniently. This being understood we will henceforth write Q(A) whenever 
the dependence of Q(A, /i) on yu is not of interest. 



3.2.2 Properties 

Theorem 1. — Properties of T- and Q-operators — 

(A) Analyticity 

The operator A^T(A) is a polynomial in A^ of degre^ N := N -f — 1 while the operator 
Q(A) is a polynomial in A of maximal degree 2/N. In the case N odd the operators Q2m '■ = 
limA->oo A^^™Q(A) and Qo := Q(0) are invertible operators and the normalization of the Q- 
operator can be fixed by Q2m = id. 

^Here, we use the notation cn = 1 for even N, cn = otherwise. 
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(B) Baxter equation 

The operators T(A) and Q(A) are related by the Baxter equation 

T(A)Q(A) = aN(A)Q(g-iA) + d^{X)Q{q\) , (3.12) 
with coefficient functions 



N 



(3.13) 



r=l 
N 

r=l 

(C) COMMUTATIVITY 

[Q(A),Q(/.)] = 0, 

VA,u. (3.14) 

[T(A),Q(/.)] = 0, 

(S) Self-adjointness 

Under the assumption C,r and Kr ^^al or imaginary numbers, the following holds: 

(T(A))t = T(A*), (Q(A))t = Q(A*). (3.15) 

For the reader's convenience we have included a self-contained proof in Appendix |Bl 

It follows from these properties that T(A) and Q(/i) can be diagonalized simultaneously for all 
A, /i. The eigenvalues Q{X) of Q(A) must satisfy 

t(A)g(A) = aN(A)g(g-iA) + dN(A)g(gA) . (3.16) 

It follows from the property (A) of Q(A) that any eigenvalue g(A) must be a polynomial of 
order 2ZN normalized by the condition Q2m = 1- Such a polynomial is fully characterized by 
its zeros Ai, . . . , A2iN, 

2ZN 

Q(A) = n(^-^'^)- (3.17) 

k=l 

It follows from the Baxter equation (13.161) that the zeros must satisfy the Bethe equations 

^i^r) ^ T-r A^ - Xrq H 1 

d(A.) l\ Xs - Xr/q ■ ^ ^ 

What is not clear at this stage is if for each solution of the Bethe equations (13.181) there indeed 
exists an eigenstate of T(A) and Q(/i). In order to show that this is the case we need a method 
to construct eigenstates from solutions to (13.181) . The Separation of Variables method will give 
us such a construction, replacing the algebraic Bethe ansatz in the cases we consider. 
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3.3 Integrability 



In order to recover the light-cone dynamics discussed in subsection l2.4l let us temporarily return 
to the homogeneous case where ^„ = 1 and k„ = k for n = 1, . . . , N. Let us note that 
the operators Y(A) simplify when A is sent to or oo. Multiplying by suitable normalization 
factors one find the unitary operators 

Yo = 7o^Y(0) and = lim 7^ /x-''^Y(/i) , 

where 70 = n!=i(l ~ ^'^^) ^^'^ 7oo = (— l)'^/' nr=i(l ~ l"^^^^)- The operators Yq and Yqo have 
the simple matrix elements 

N 



(z|Yo|z') = l[q 



Y^|z') = l[q 



n=l 
N 




if <( ^ } (3.19) 



n=l 

and 



Q+(A) = Y(A)-YL, Q-(A) = (Y(A)-Y5)-' (3.20) 
Integrability follows immediately from the following observation: 



U = U+ ■ U-, U+ = Q+(l/«e), U- = Q-(«:/e), 



(3.21) 



where = Y[l=ii^ ^ q'^'^~^Y^ / {k^ — q'^^'^^Y^ . The proof can be found in Appendix B. 
It is very important to remark that there is of course no problem to construct time evolution 
operators in the inhomogeneous cases by specializing the spectral parameter of the Q-operator 
in a suitable way. We are just not able to represent the time evolution as simple as in (I2.14I ). 
One will still have a lattice approximation to the time evolution in the continuum field theory as 
long as the inhomogeneity parameters are scaled to unity in the continuum limit. 



4. Separation of variables I — Statement of results 

The Separation of Variables (SOV) method of Sklyanin llSk85ll - [lSk95ll as developed for lattice 
Sine-Gordon model in this section will allow us to take an important step towards the simulta- 
neous diagonalization of the T- and Q-operators. 

The separation of variables method is based on the observation that the spectral problem for 
T(A) simplifies considerably if one works in an auxiliary representation where the commutative 
family B(A) of operators introduced in (13.31) is diagonal. In the following subsection we will 



16 



discuss a family of representations of the Yang-Baxter algebra (13.41) that has this property. We 
will refer to this class of representations as the SOV-representations. We will subsequently 
show that our original representation introduced in (|3.1|) . (|3.2I) is indeed equivalent to a certain 
SOV-representation. 



4.1 The SOV-representation 

The operators representing (13.41) in the SOV-representation relevant for the case of a lattice with 
N sites will be denoted as 

We will now describe the representation of the algebra (13.41) in which Bn(A) acts diagonally. 



4.1.1 The spectrum of 

By definition, we require that Bn(A) is represented by a diagonal matrix. In order to parame- 
terize the eigenvalues, let us fix a tuple C, = (C^, . . . , C,^) of complex numbers such that C'p ^ 
for a ^ h. The vector space C^^ underlying the SOV-representation will be identified with the 
space of functions ^(r/) defined for rj taken from the discrete set 

Bn = { (g'^Ci, • • • , g'"CN) ; (fci, . . . , ^n) e } . (4.2) 

The SOV-representation is characterized by the property that B(A) acts on the functions '^{rf), 
T] = (r^i, . . . , ?7n) G Bn as a multiplication operator, 

{X/Va - Va/X) ; (4.3) 
n=l a=l 

where [N] = N — cn- We see that r]i, . . . , r]^^ represent the zeros of 6,, (A). In the case of even 
N it turns out that we need a supplementary variable ?7n in order to be able to parameterize the 
spectrum of B(A). 



Bn(a) ^iv) = MX) ^iv) , MX) ^ n - 



4.1.2 Representation of the remaining operators 

Given that Bn(A) is represented as in (14.31 ). it can be shown ||Sk85ll - IISk95lFI that the represen- 
tation of the remaining operators An (A), Cn(A) Dn(A) is to a large extend determined by the 



See liBT091 for the case of the Sinh-Gordon model which is very similar to the case at hand. 
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algebra (13.41) . First note (see e.g. IIBT091 Appendix C.2] for a proof) that the so-called quantum 
determinant 

detq(M(A)) = A(A)D(g-^A) - B(A)C(g"^A) (4.4) 

generates central elements of the algebra (13.41) . In the representation defined by (13. Ik (13.21 ) we 
find that A^'^detq(M(A)) is a polynomial in A^ of order 2N. We therefore require that 

AN(A)DN(g-U) - BN(A)CN(g-^A) = An(A) ■ id , (4.5) 

with A^'^An(A) being a polynomial in A^ of order 2N. 

The algebra (|3.4I) furthermore implies that An (A) and Dn(A) can be represented in the form 



An(A) = e^br,{X) 



Dn(A) = ej^briiX) 



Va 



Va 



t; 



At- - 



N 



[N] 

a=l b=ia 
[N] 



a=l by^a 



X/Vb - Vb/ A 

Va/Vb- Vb/Va 

A/ r]b - Vb/ A 

Va/Vb - Vb/Va 



(4.6) 
(4.7) 



where are the operators defined by 



The expressions (14.61) and (14.71) contain complex- valued coefficients r^^, r]^, a]<f{r]r) and dt^{r]r 
The coefficients ajsi{r],r) and d^(r]r) are restricted by the condition 



,N, 



(4.8) 



as follows from the consistency of (14.51) . (14.31) . (14.61) and (|4.7I) . This leaves some freedom in the 
choice of a]<i{r]r) and d]<i{r]r) that will be further discussed later. 

The operator Cn(A) is finally univocall)@ defined such that the quantum determinant condition 
(1431) is satisfied. 



4.1.3 Central elements 

For the representations in question, the algebra (|3.4I) has a large center. For its description let 
us, following ||Ta91|| . define the average value O of the elements of the monodromy matrix 

Ms°^(A) as 

p 

0(A) = nO(g'A), A = A^ (4.9) 

k=l 

where can be An, Bn, Cn or Dn. 

^Note that the operator Bn (A) is invertible except for A which coincides with a zero of Bn, so in general Cn (A) 
is defined by (4.5) just inverting Bn(A). This is enough to fix in an unique way the operator Cn being it a Laurent 
polynomial of degree [N] in A. 
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Proposition 1. The average values ^n(A), B^{K), Cn(A), V^{K) of the monodromy matrix 
M(A) elements are central elements. 

The Proposition is proven in [|Ta91| . see Subsection 15.21 for an alternative proof. The average 
values are of course unchanged by similarity transformations. They therefore represent pa- 
rameters of the representation. Let us briefly discuss how these parameters are related to the 
parameters of the SOV-representation introduced above. 

First, let us note that B^{k) is easily found from (|4.3I) to be given by the formula 

N 7 — P 

^" "^(A/Z, - VA) , (4.10) 



n=l " a=l -"-a — "'a 

The values A^^iZ^) and 'D^(Zr) are related to the coefficients a^{q^rir) and d^{q^rir) by 

p p 
= n^N(gV), V^{Z,) = Wd^{q''r]r) . (4.11) 

k=l k=l 

Note that the condition (14.81) leaves some remaining arbitrariness in the choice of the coefficients 
a^{'q), d]ss{ri). The gauge transformations 

N 
r=l 

induce a change of coefficients 

ONl^r) = aniVr) j,^^^ , d^iVr) = "N(r/r) ^ , (4.13) 

but clearly leave A^{Zr) and V^(Zr) unchanged. The data A^{Zr) and V^{Zr) therefore 
characterize gauge-equivalence classes of representations for An (A) and Dn(A) in the form 
(1431) . 



4.2 Existence of SOV-representation for the lattice Sine-Gordon model 

We are looking for an invertible transformation W^°^ that maps the lattice Sine-Gordon model 
defined in the previous sections to a SOV-representation, 

(W'°^)-^ ■ M'°^(A) ■ W"°^ = M(A) . (4.14) 

Constructing M^°^(A) is of course equivalent to the construction of a basis for H consisting of 
eigenvectors (77 1 of B(A), 

(r;|B(A) = v'-b,{X){v\. (4.15) 
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The transformation W^°^ is then described in terms of ( 77 1 2; ) as 

(W^°»(7^) = (^I^)^W- (4-16) 

2G{Sp)N 

The existence of an eigenbasis for B(A) is not trivial since B(A) is not a normal operator. It 
turns out that such a similarity transformation exists for generic values of the parameters u, v, ^ 
and K. 

Theorem 2. - Existence of SOV-representation for the lattice Sine-Gordon model - 

For generic values of the parameters u, v, C, and k there exists an invertible operator W^°^ : 
n which satisfies KM . 

The proof is given in the following Section[5l It follows from (14.61) . (14.71) that the wave-functions 
\E'(?7) = ( ?7 1 1 ) of eigenstates 1 1 ) must satisfy the discrete Baxter equations 

t(r70*(^) = a(r/„)T;*(r7) + rf(r/„)T+^(r/) , (4.17) 

where n = 1, . . . , N. Equation (14.171) represents a system of linear equations for the p^ 
different components "^{r]) of the vector \E'. It may be written in the form Dt-'^ = 0, where Df 
is a p'^xp^ -matrix that depends on t = t{\). The condition for existence of solutions detD^ = 
is a polynomial equation of order p^ on t{X) . We therefore expect to find p^ different solutions, 
just enough to get a basis for V.. 

We will return to the analysis of the spectral problem of T(A) in Section[6l Let us now describe 
more precisely the set of values of the parameters for which a SOV-representation exists. 



4.3 Calculation of the average values 

Necessary condition for the existence of W^°^ is of course the equality 

M{A) = M'°''(A), (4.18) 

of the matrices formed out of the average values of M(A) and M^°^(A), respectively. It turns 
out that Ai (A) can be calculated recursively from the average values of the elements of the Lax 
matrices L^{\), which are explicitly given by 

^ (jy. ]^( ^'u^iKVn + V-') K„(Av;/x„ - x„/KA)\ 

where we have used the notations Kn = k,^, X„ = ^p, Un = u"^ and Vn = v^. Indeed, we have: 
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Proposition 2. We have 

a^n(a) = /:n(a) /:n-i(a) . . . /:i(a) . (4.20) 

This has been proven in [|Ta91| . see Subsection 15 .21 for an alternative proof. 

The equality (14.181 ) defines the mapping between the parameters u,v, k and ^ of the representa- 
tion defined in Subsection [3T| and the parameters of the SOV-representation. Formula (|4.20l) in 
particular allows us to calculate B{A) in terms of u, v, k and ^. Equation (|4.10l) then defines the 
numbers Za = r]^ uniquely up to permutations of a = 1 , . . . , [N] . 

Existence of a SOV-representation in particular requires that Za 7^ ^6 for all a ^ b, a,b = 
1, . . . , [N]. It can be shown (see Subsection 15 .3l below) that the subspace of the space of param- 
eters u, V, K and ^ for which this is not the case has codimension at least one. Sufficient for 
the existence of a SOV-representation is the condition that the representations T^m exist for all 
M = 1,...,N-1. 



5. Separation of variables II — Proofs 



We are now proving Theorem |2] by constructing a set of linearly independent vectors ( 77 1 
which are eigenvectors of B(A) with distinct eigenvalues. This will be equivalent to a recursive 
construction of the matrix of elements {ri\z) and so of the invertible operator W^°^ : "H — )■ 
^sov relation (I4l6l) . 



5.1 Construction of an eigenbasis for B(A) 

We will construct the eigenstates (?7 1 of B(A) = Bn(A) recursively by induction on N. The 
corresponding eigenvalues B{\) are parameterized by the tuple rj = (r7a)a=i n 

N [N 

B{X) = r/^- 6,(A) , 6,(A) = J] " H (^/^« ' ^'^Z^) ' (^.1) 



2 

n=l a=l 

We remind that is zero for N odd and 1 for N even. 

In the case N = 1 we may simply take ( r/i | = ( f | , where ( f | is an eigenstate of the operator 
vi with eigenvalue v. It is useful to note that the inhomogeneity parameter determines the subset 
of C on which the variable rji lives, r]i G ^iSp. 

Now assume we have constructed the eigenstates ( x I of Bm(A) for any M < N. The eigenstates 
( ?7 1, ?7 = (r^N, . . . , r]i), of Bn(A) may then be constructed in the following form 

^^1 = ^n(^Ixz;Xi)(XzI®(XiI , (5.2) 

Xi Xz 
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where ( Xz I and ( | are eigenstates of Bm(A) and Bn_m(A) with eigenvalues parameterized as 
in dSU) by the tuples Xz = {Xza)a=i,...,M and Xi = (Xia)a=i,...,N-M' respectively. It suffices to 
consider the cases where N — M is odd. 



It follows from the formula 

Bn(A) = Am (A) ® B (A) + Bm(A) g) D 

N— M (A) 

= A,m(A)B 

1 N— M 

(A) + B,m(A)D 

1 N— M (A) 

that the matrix elements K,^( rj \ x^] Xi ) have to satisfy the relations 
A2m(A)B 

1 N— M 

(A) + B,m(A)D 

1 N— M V I Xi, Xi ) 



(5.3) 



N [N] 



(5.4) 



= n f n (^/^'^ - ^'^z^) ^ I X.; X. ) , 

n=l a=l 

where we used the notation 0* for the transpose of an operator 0. 
Let us assume that 

X.aq"' i A„ x^bq''^' i A, and x.aq""' ^ X^bq"', 

where hi E {1, ...,p}, a E {1, ...,N — M} , b E {I, ...,M} and A, is the set of zeros of the 
quantum determinant on the subchain i, with i = 1,2. Under these assumption^ the previous 
equations yield recursion relations for the dependence of the kernels in the variables Xia and 
Xib simply by setting A = Xia and A = Xib- Indeed for A = Xia the first term on the left of (|5.4I) 
vanishes leading to 



(5.5) 



N-M . [M] 



'^iaK^{v\Xz,xJ d^{q \J x7 n —HixJXzb-xJxJ 



n=l a=l 
[Nl 



/n 



(5.6) 



ixJVb-Vb/x la J 



6=1 



= K^iv\Xz;Xi) vt 

while for A = Xm one finds similarly 

'^laKAv\Xz;Xj a^iq^^za) n~ YliXza/Xib-XjXz 



M . N-M 



1 I 1 

n=l 6=1 



(5.7) 



iXia/Vb-Vb/X 1 



6=1 



If M is even we find the recursion relation determining the dependence on Xim by sending 
A — )■ 00 in (15.41) . leading to 

^ M-l ^ N-M ^ N ^ 

tL^n(^Ix.;xJ — n — n — = ^^iv\xz;v.) U- (^.8) 

-I Xib 



XzA ^—i Xza 



= 1^^> 



The subspace within the space of parameters where these conditions are not satisfied has codimension at least 



one. 



N-M .p [M] 
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The recursion relations (15.61) . (I5.7|) have solutions compatible with the requirement of cyclicity, 
{Tj^^y = 1 and (T^^)^ = 1 for all values of a, provided that the algebraic equations 

[Nl 

'ixUv! - </xla) , 

(5.9) 

where £'i(Xia) = JJ^^^^^'^i")' 

k=l 

and 

M .p N-M N 

X{]jX{{xlJxl,'xlJxla) = ivTYU^xUvl-vUxla), 



n=l b=l 6=1 



n=l b=l b=l 

P 

where A^{xia) = Yla^{q''xia) , 

k=l 

are satisfied. If M is even the recursion relation (15.81) yields the additional relation 

M-l N-M N 



(5.10) 



^2A a=l ^^'^ 6=1 Vb 



(5.11) 



We will show in the next subsection that the equations (|5.9I )- (I5.1 II) completely determine 7]^ in 
terms of xL' xL- 

By using (14.101 ) and (16.91 ) it is easy to see that the conditions (15.91 ) and (15.101 ) are nothing but the 
equations 

B^{A) = ^m(A)Sn-m(A) + i3M(A)2^N-M(A), (5.12) 

evaluated at A = Xia A = xla^ respectively. The relation (15.1 II) follows from (15.121) in the 
limit A — )■ oo. The relations (|5.12l) are implied by (14.201) . We conclude that our construction 
of B(A)-eigenstates will work if the representations TZj^, TZm and 7^n-m are all non-degenerate. 
Theorem 2 follows by induction. 



5.2 On average value formulae 

Proposition 3. The average values of the Yang-Baxter generators are central elements which 
satisfy the following recursive equations: 

Bn{A) = ^M(A)SN-M(A) + i3M(A)PN-M(A), (5.13) 

Cn(A) = Pm(A)Cn-m(A)+Cm(A)^n-m(A), (5.14) 

^n(A) = ^M(A)^N-M(A) + i3M(A)CN-M(A), (5.15) 

I?n(A) = I?M(A)I)N-M(A)+CM(A)i3N-M(A), (5.16) 

where N — M or M odd. 
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Proof. In the previous subsection we have proven the existence of SOV-representations, i.e. 
the diagonalizability of the B -operator. First of all let us point out that A(A), B(A), C(A) and 
D(A) are one parameter families of commuting operators. This implies that the corresponding 
average values are functions of A = A^. 

The fact that B^{K) is central trivially follows from the fact that Bn(A) is diagonal in the 
SOV-representation, while for the operators A and D we have that for N odd, ^n(A)A^^^ and 
X'n(A)A^^^ are polynomials in A^ of degree N — 1. It follows that the special values given by 
(14.111 ) characterize them completely, 



where AN(^a) and D^{Za) are the average values of the coefficients of the SOV-representation. 
In the case of N even we have just to add the asymptotic property of ^n(A) and 'D^{K) dis- 
cussed in appendix O to complete the statement. Finally, the fact that Cn(A) is central follows 
by its diagonalizability in the cyclic representations. 

Now the above recursive formulae (|5.13H5.16l) are a simple consequence of the centrality of the 
average values of the monodromy matrix elements. Let us consider only the case of the average 
value of An (A). We have the expansion: 



in terms of the entries of the monodromy matrix of the subchains i and 2 with (N — M) -sites 
and M-sites, respectively. It follows directly from definition (14.91 ) of the average value together 
with (15.181) that ^n(A) can be represented in the form 



where An (A) is a sum over monomials which contain at least one and at most p — 2 factors 
of A2 M(Ag™). As before, we may work in a representation where the B2 m{M") are diagonal, 
spanned by the states {xi \ introduced in the previous subsection. As the factors A^ uiM^) 
contained in An (A) produce states with modified eigenvalue of B^ m{M^)^ none of the states 
produced by acting with An (A) on ( I can be proportional to {Xi\- This would be in contra- 
diction to the fact that ^n(A) is central unless An(A) = 0. □ 




(5.17) 




rtN 

(A) = A,m(A)A, (A) + B,m(A)Q N-M (A), 



(5.18) 



^n(A) = Am(A)A N-M N-M (A) + An(A) 



(5.19) 
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5.3 Non-degeneracy condition 

Proposition 4. The condition Zr = Zg for certain r ^ s with r,s G {1, [N]} defines a 
subspace in the space of the parameters {ki, kn, ^n} G of codimension at least 
one. 



Proof. The parameters Zr are related to the expectation value B^{A) by means of the equation 

\A/Z, - ZJA) . (5.20) 



n=l a=l 

It follows from (14.201) and (14.191 ) that B^{K) is a Laurent polynomial in X„ that depends poly- 
nomially on each of the parameters Kn- Equation (15.201) defines the tuple Z = (Z^, . . . , Z^^^) 
uniquely up to permutations of Z^, . . . , Z^^^ as function of the parameters X = (Xi, . . . , Xn) 
and K = {Ki, . . . , K-!<s). We are going to show tha|zl 

JiX;K) = det(^) ^0. (5.21) 

\0^s/ r,s=l,...,[N] 

The functional dependence§ of the Z^, . . . , Z^-^-^ w.r.t. the parameters K implies that it is suffi- 
cient to show that J{X; K) ^ for special values of K in order to prove that J{X; K) ^ 
except for values of K within a subset of of dimension less than N. 

Let us choose Kn = for n = 1, [N], then the average values (14.191 ) of the Lax operators 
simplify to 

£-(A)=f , ° , '^A"-^-»/M. (5.22) 
Inserting this into (14.201) yields 

[N 

Bn(A) = {Kl, + 1)^^* HiA/X^ - XJA) . (5.23) 

n=l 

The fact that J{X; K) follows for the case under consideration easily from (15.231) . 

Whenever J{X; K) ^ 0, we have invertibility of the mapping Z = Z{Xi, . . . , X[n]). The claim 
follows from this observation. □ 



^It should be noted that for even N it is indeed sufficient to consider the dependence w.rt. Xi, . . . , Xn-i 
'^Let o-lf^l (Z) be the de gree n elementary symmetric polynomial in the variables Z, then al^^ (Z)/ CTjJ^j (Z) are 
Laurent polynomials of degree 1 in all the parameters X and K 



6. The spectrum — odd number of sites 
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Let us now return to the analysis of the spectrum of the model. For simplicity we will consider 
here the case of odd N, while we will discuss the case of even N in the next section. The 
existence of the SOV-representation allows one to reformulate the spectral problem for T(A) as 
the problem to find all solutions of the discrete Baxter equations (|4.17|) . This equation may be 
written in the form 

Vr ^(r^) = , Vr= a{r]r)T; + d{r]r)Tt - , (6.1) 
where r = 1, . . . , N. Previous experience with the SOV method suggests to consider the ansatz 

N 
r=l 

where Qt{^) is the eigenvalue of the corresponding Q-operator which satisfies the- functional 
Baxter equations 

t{X)Qt{\) = M>^)Qtiq-'X) + d^{X)Qt{qX) . (6.3) 

This approach will turn out to work, but in a way that is more subtle than in previously analyzed 
cases. 

6.1 States from solutions of the Baxter equation 

First, in the present case it is not immediately clear if the functional Baxter equation (|6.3I) and 
the discrete Baxter equation (|6.1I) are compatible. The question is if one can always assume that 
the coefficients a(r]j.) and d{r]r) in (16.11) coincide with the coefficients a^(r]r), d^ir],.) appearing 
in the functional equation (16.31) satisfied by the Q-operator. The key point to observe is contained 
in the following Lemma. 

Lemma 1. Let Aj^(A) and Dj^(A) be the average values of the coefficients aN(A) and d^{X) of 
the Baxter equation f |(5.iD . 

V P 

An(A) = n^N(g'A), Dn(A) = ndN(g'A). (6.4) 

fc=i fc=i 

We then have 



An (A) = ^N(A)-i3N(A), Dn(A) 



= ^n(A) + Sn(A) 



(6.5) 
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Proof. The claim is checked for N = 1 by straightforward computation. Let us assume now 
that the statement holds for N — 1 and let us show it for N. The average values An (A) and Dn(A) 
satisfy by definition the factorization: 

An(A) = Af )(A)A{f_-/'-'^)(A), Dn(A) = of )(A)D{f_-^'-'^)(A), (6.6) 

where the upper indices are referred to the quantum sites involved while the lower indices to 
the total number of sites. We can use now the induction hypothesis to get the result: 

An(A) = (4^) - BT\mA^:^-'\k) - <_V'-'')(A)) = ^n(A) - Sn(A), (6.7) 
Dn(A) = (4^) + BT\k))US-i'-''\^) + <-V'-''^(A)) = ^n(A) + Sn(A), (6.8) 

where in the last formulae we have used (14.201) together with the fact that .4n(A) = V^[Is) and 
Bn(A) = Cn(A) for u„ = 1, = 1, n = 1, . . . , N. □ 

The Lemma implies in particular 

k^{Zr) = A^{Zr) , Dn(Z,,) = V^{Z,) , (6.9) 

for all r = 1, . . . , N. We may therefore always find a gauge transformation (14.131) such that the 
coefficients a^{rij) and dj^{r]r) in (16.11) become equal to 

%('7r) = aN('7r) , dj^{r]r) = dj^{r]r) , (6.10) 

respectively. So from now on we will denote also the coefficients in (|4.17l) with a and d omitting 
the index N unless necessary. The ansatz (|6.2I) therefore indeed yields an eigenstate of T(A) for 
each solution Qt{\) of the functional Baxter equation (|3.16l) . We are going to show that all 
eigenstates can be obtained in this way. 

6.2 Non-degeneracy of T( A) -eigenvalues 

In order to analyze the equations (|6.1I) . let us note that the matrix representation of the operator 
Vr defined in (16.11) is block diagonal with blocks labeled by n = 1, . . . , N. Let '^niv) ^ be 
the vector with components 

"^nAv) = '^iVl, • • • , Vn-1, Cnq^, Vn+1, ■ ■ ■ , Vn) ■ 

Equation (|6.1I) is then equivalent to the set of linear equations 



L)W.^^(ry) = 0, r = l,...,N. 



(6.11) 
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where D^"^^ is the p x p-matrix 

/ t{Cr) -d(Cr.) 

-a(gCr) t{qCr) -d(gCr 




-a(Cr) \ 








-a(g2'C) 



^21, 



(6.12) 



... 

\-d{q''Cr) ... -a(g2'C) t{q''Cr) J 

The equation (16.111 ) can have solutions only if det{D^'^^) = 0. The determinant det{D^'^^) is a 
polynomial of degree p in each of the N coefficients of the polynomial t(A). 

Proposition 5. G/ven that det (D^'') ) = 0, the dimension of the space of solutions to the equation 
A6.11\i for any r = 1, ... ,1^ is one for generic values of the parameters ^ and n. 

Proof. Let us decompose the p x p matrix D^'^'> into the block form 



(6.13) 



where the submatrix E^'^'^ is a (p — 1) x (p — 1) matrix, t;^''^ and w^^^ are column and row vectors 
withp — 1 components, respectively. We assume that det(D^'"^) = 0, so existence of a solution 
to \1/ = is ensured. It is easy to see that the equation Z}('")\[^ = has a unique solution 
provided that det{E^''^) ^ 0. 

It remains to show that det{E^'^'^) ^ holds for generic values of the parameters ^ and k. To 
this aim let us observe that the coefficients a.(q''(j.) and d(q''(j.) appearing in (16.1 II ) depend 
analytically on the parameters k. If det{E^^^) = is not identically zero, it can therefore only 
vanish at isolated points. It therefore suffices to prove the statement in a neighborhood of the 
values for the parameters k which are such that 



a(C) = 

Such values of k and ^ exist: Setting 



d{q-\r 



0. 



±i forn = 1, . . . , N, one finds that 



n(A/^n-X/A) , 



(6.14) 



(6.15) 



n=l 



which vanishes for A = g2^„. We may therefore choose§ (n = Q^^n- We then find (16.141 ) from 
(I3l3]). (lOOl). 

''Note that this choice impHes that v„ e (-Ij^'/^gr^/^Sp. 
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Given that (16.141) holds, it is easy to see that det(-E''^^^) 7^ . Indeed, the submatrix E^i \ is lower 
triangular if (16.141) is valid, and it has — d(g'^^^.). A; = 0, ... ,p — 2 as its diagonal elements. It 
follows that det{E^^^) = Y[k=o ^(^'^O which is always nonzero if (16.141) is satisfied. □ 

The previous results admit the following reformulation which is central for the classification 
and construction of the spectrum of T(A): 

Theorem 3. For generic values of the parameters k and ^ the spectrum o/T(A) is simple and 
all the wave-functions "^tiv) '^'^^ represented in the factorized form A6.2\l with Qt being the 
eigenvalue of the Q-operator on the eigenstate \ t). 

The eigenvectors \ t) o/T(A) are in one-to-one correspondence with the polynomials QtW of 
order 2/N, with Qt{0) 7^ 0, which satisfy the Baxter equation ( 13.7(51) with t{X) being an even 
Laurent polynomial in A of degree N — 1. 

Proof. Proposition [5] implies that the spectrum of T(A) is simple. Let It) be an eigenstate 
of T(A). Self-adjointness and mutual commutativity of T(A) and Q(/u) imply that 1 1) is also 
eigenstate of Q(A). Let Qt{\) be the Q-eigenvalue on \ t). The polynomial Qt{X) is related to 
t(A) by the Baxter equation (|3.16l) which specialized to the values \ = r]r yields the equations 
(|6.1 II) . It follows that there must exist nonzero numbers u^. such that 

QtiCq") = z/.^.,fc(Ci,...,CN). (6.16) 

This implies that the wave-functions \&(r/) can be represented in the form (|6.2I) with Qt being 
the eigenvalue of the Q-operator on the eigenstate \t). □ 

Remark 1. It may be worth noting that the equivalence with the Fateev-Zamolodchikov model 
does not hold for odd number of lattice sites. The spectrum of the two models is qualitatively 
different, being doubly degenerate in the Fateev-Zamolodchikov model but simple in the lattice 
Sine-Gordon model, as illustrated in Appendix iDl 



6.3 Completeness of the Bethe ansatz 

Assume we are given a solution (Ai, . . . , \2in) of the Bethe equations (13.181 ). Let us construct 
the polynomial Q{\) via equation (13.171) . Define 

i^X) :- ^(A)Q(g-^A) + d(A)g(gA) _ ^^^^^ 
Q(A) 

t(A) is nonsingular for A = A^, /c = 1, . . . , M thanks to the Bethe equations (13.181) . The pairs 
iQii]r),t{r]r)) satisfy the discrete Baxter equation by construction. Inserting this solution into 
(16.21) produces an eigenstate 1 1 ) of the transfer matrix T(A) within the SOV-representation. 
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Conversely, let 1 1 ) be an eigenvector of T(A) with eigenvalue t(A). Let Q't{X) be the eigenvalue 
of Q(A) on \ t). Thanks to the properties of Q(A) listed in Theorem [Done may factorize QJ(A) 
in the form (13.171) . The tuple of zeros (A'^^, . . . , \'2m) of Q't{\) must satisfy the Bethe equations 
(13.181) as follows from the Baxter equation (|3.12l) satisfied by Q(A). Inserting Q't{'qr) into (|6.2I) 
produces an eigenstate 1 1' ) that must be proportional to 1 1 ) due to the simplicity of the spectrum 
ofT(A). 

It follows that there is a one-to-one correspondence between the solutions to (|3.18|) and the 
eigenstates of the transfer matrix {Completeness of the Bethe ansatz). 



7. The spectrum — even number of sites 



We will now generalize these results to the case of a chain with even number of sites. It 
turns out that the spectrum of T(A) is degenerate in this case, but the degeneracy is resolved by 
introducing an operator 9 which commutes both with T(A) and Q(A). The joint spectrum of 
T(A), Q(A) and 6 is found to be simple. 



7.1 TheG-charge 



In the case of a lattice with N even quantum sites, we can introduce the operator: 

N 

e = nv(r^^""- (7.1) 

n=l 

Proposition 6. O commutes with the transfer matrix and satisfies the following commutation 
relations with the entries of the monodromy matrix: 

eC(A) = gC(A)a [A(A),e] = 0, (7.2) 

B(A)0 = g0B(A), [D(A),e]=0. (7.3) 



Proof. The claim can be easily verified explicitly for N = 2. The proof for the case of general 
even N = 2M follows by induction. Indeed, 



A22M Ai2(N-M) + B22M Ci2(N-M) A22M Bi2(N-M) + B22M Di2(N-M) 
C22M Ai2(N-M) + D22M Ci2(N-M) C22M Bi2(N-M) + D22M Di2(N-M) 



M22M Mi2(N-M) = 

which easily allows one to deduce that the claim holds if it holds for all M < N. □ 
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7.2 T-9-spectruin simplicity 

Lemma 2. Let k G {— /, .., /} and \tk) be a simultaneous eigenstate of the transfermatrix T(A) 
and of the Q-charge with eigenvalues t\k\ (A) and q^, respectively, then X^t\k\ (A) is a polynomial 
in A^ of degree N which is a solution of the system of equations: 

det(D(")) = Vr G {1,..., [N]}, (7.4) 

where the p x p matrices D^*") are defined in ^6.12\l , with asymptotics oft\k\ (A) given by: 

^ lim A^^t|,|(A) = fn [q' + q-') . (7.5) 

logA-i.±oo V 1 ^ / 

Proof. The fact that the generic eigenvalue of the transfer matrix has to satisfy the system 
(17.41) has been discussed in Section[6l so we have just to verify the asymptotics (17.51) for the T- 
eigenvalue t^^i i^)- This follows by the assumption that \tk) is an eigenstate of with eigenvalue 
g*^, and by formulae 

lim A^^T(A) = I n (e + , (7.6) 

logA->±oo \ 1 ^ / 

derived in appendix O □ 
The previous Lemma implies in particular the following: 

Theorem 4. For generic values of the parameters k and ^ the simultaneous spectrum ofT and 
6 operators is simple and the generic eigenstate of the T-Q-eigenbasis has a wave-function 
of the form 

N-l 

'^{v)=VN'll'^\k\{Va), (7.7) 

a=l 

where, for any r G {1,...,N — 1}, the vector (V^ifcilCr), V'|fc|(CrQ'), V'|A:|(Cr'?^0) unique 
(up to normalization) solution of the linear equations d<5.iiD corresponding to t|fc|(A). 

Proof. Let us use the SOV-construction of T-eigenstates and let us observe that an analog of 
Proposition [5] also holdil*] for even N. This implies that the wave-function '^{rf) can be repre- 
sented in the form 

N-l 

a=l 

Finally, using that \tk) is eigenstate of with eigenvalue q^ we get ft^ijl^) oc rj^^ . □ 



10' 



The proof given previously holds for both the cases N even and odd just changing N into [N] everywhere. 
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Thanks to the explicit construction of the simultaneous T-0 eigenstates given in (|7.7I) . we have 
that the eigenstates of T(A) with 6-charge eigenvalue 1 are simple, while all the others are dou- 
bly degenerate with eigenspaces generated by a pair of T-eigenstates with 6-charge eigenvalues 

7.3 Q-operator and Bethe ansatz 

Let us point out some peculiarity of the Q-operator in the case of even chain. In order to see 
this, we need the following Lemma which is of interest in its own right. 

Lemma 3. For a given t{\), there is at most one polynomial of degree 2/N which satisfies the 
Baxter equation A3.16\) . 

Proof. Let us define the q-Wronskian: 

W{\) = Qii\)Q2{q-'\) - Q2{\)Qiiq-'\) . (7.9) 

written in terms of two solutions Qi{X) and Q2(A) of the Baxter equation; then W{\) satisfies 
the equation 

a(A)Py(A) = d(A) T+1^(A) . (7.10) 
Note now that Lemma [T] implies: 

21 21 

J]a(Ag'=) ^ VA^Bn, (7.11) 

k=0 k=0 

SO for any A ^ Bn the only solution consistent with cyclicity (T"'")^ = 1 is W{\) = 0. It is then 
easy to see that this implies that Qi{\) = Q2i^)- CH 

Now we can prove the following: 

Proposition 7. The Q-operators commute with the Q-charge and \t±\k\) are Q-eigenstates with 
common eigenvalue Q\k\ (A) of degree 2/N — k{a^p ± 1) in A and a zero of order k{a^p ± 1) at 
A = 0, where and a'^ are non-negative integers, while Oq and are positive integers. 

Proof. The commutativity of T and Q-operators implies that the T-eigenspace C{\t±^k\)) cor- 
responding to the eigenvalue t|fc|(A) is invariant under the action of Q and so for /c = any 
T-eigenstate |to) is directly a Q-eigenstate. Let us observe that the self-adjointness of Q implies 
that in the two-dimensional T-eigenspace C{\t±\k\)) with A; 7^ we can always take two linear 
combinations of the states |t|fe|) and |t-|A;|) which are Q-eigenstates. Now thanks to the Lemma 
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[3] for fixed T-eigenvalue t|fe|(A) the corresponding Q-eigenvalue Q\k\i^) is unique which im- 
plies that |t±|fc|) are themselves Q-eigenstates. The commutativity of the Q-operator with the 
9-charge follows by observing that the \t±\k\) define a basis. 

Let us complete the proof showing that the conditions on the polynomial Q\k\i^) stated in the 
Proposition are simple consequences of the fact that |t±|fc|) are eigenstates of the 6-charge with 
eigenvalues Indeed, the compatibility of the asymptotics conditions (17.51) with the TQ 

Baxter equation implies 

lto«H^ = ,±l'-l ,toM^.,-(N±™, (7.12) 

which are equivalent to the conditions on the polynomial Q\k\ (A) stated in the Proposition. □ 

Note that the uniqueness of the Q-eigenvalue Q\k\{^) corresponding to a given T-eigenvalue 
t|fe|(A) implies that each vector {il'\k\iCr),'ip\k\iCrC[), ■■■ii>\k\{C,r<f^)) appearing in (|7.7I) must be 
proportional to the vector {Q\k\{C,r), Q\k\iCrC[), Q\k\{C,r(f^)) so that the previous results admit 
the following reformulation: 

Theorem 5. The pairs of eigenvectors \t\k\) and \t^\k\) ofT{X) are in one-to-one correspon- 
dence with the polynomials Q\k\{X) of maximal order 2/N which have the asymptotics A7.12\) 
and satisfy the Baxter equation d3.i(5D with t\k\{\) being an even Laurent polynomial in A of 
degree N. 

As in the case of N odd this reformulation allows the classification and construction of the 
spectrum of T(A) by the analysis of the solutions to the system of the Bethe equations. 



A. Cyclic solutions of the star-triangle relation 

It will sometimes be convenient for us to identify Zp = Z/pZ with the subset §p = {g^"; n = 
—I, ... J} of the unit circle since = 1. 



A.l Definition and elementary properties 

A.1.1 The function w\{z) 

Let us define a function : Sp — ^ C by 

^^(f-^^U-^S^UT^- " = 0,...,p-l. (A.1) 

r=l r=l 
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This function is indeed cyclic (defined on §p) since 11^=1(1 "~ = 1 — implies 

wx{q^^) = ^A(g''+') = ^a(I) , (A.2) 
The function w\{z) is the unique solution to the functional equation 

{z + \)wx{qz) = {l + Xz)wx{q-^z), (A.3) 
which is a polynomial of order / in A and which satisfies the normalization condition 

= 1 Vn e Zp. (A.4) 
The function wx{z) satisfies the inversion relation 

2r-l\{^ I \^2r-l^ 



A.1.2 The function w\{z) 

Let us also introduce the function w\{z) as the discrete Fourier transformation of w\, 

1 ' 

^a(^) = -Tz^'w^iq'') (A.6) 

wa(2;) can be characterized as the unique solution to the functional relation 

(1 - Xqz)Wx{qz) = {z - qX)wx{q~^z) , (A.7) 

which is a polynomial of order I in A and which satisfies the normalization condition Wi{q") — 
Sn,o- It may therefore be represented by the product 

r=l s=l 

It is also useful to observe that wx and wx are related by complex conjugation as follows: 

{w,x{z)r = W^z) n ^ ~ Li ■ (A.9) 

s=l 

This relation makes it easy to deduce properties of wx from those of wx- 
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A. 1.3 Further functional relations 

Let us list further functional relations satisfied by the function wx{z). 

iX + z)w^{qz) = /~^'z+^l + q\)Wg^iz), 
{1 + X)w^{qz) = q'"+'z-Hl + Xz)w^^^iz), 
{1 - qX)w^{qz) = q'^^'^ z'^z - qX)Wg^{z), 
(1 - qXz) w^iqz) = q-^'-^ z+^ (1 - A) w^,^{z). 

These relations play a key role in the derivation of the Baxter equation (|3.12l) . 



(A. 10) 



A.2 Star-triangle relation 

One of the most important properties of the function w\{x) is the star-triangle relation IIFZ82II 

Wa{x/u)Wap{x/v)wp{x/w) = W a{w / v) W ajiiu / w) W ^{v / u) , (A.U) 

see P BaOSII for an elegant proof and references to related work. We are mainly going to use the 
following consequence of (lA.l II) called the exchange relation 

^ w^{y/u)wp{y/v)w^{y/w)w^{y/x) = (A.12) 

W^, (li/v) . ^ _ _ 

= — rr^z^ ^piyl^)^c.{yh)ws{y/w)w^{y/x) . 

for a'j / (35 = 1. In order to prove (IA.12I) let us note the relation 
_ _ 1 ' 

Wa{u/z)wi/a{z/v) = - "^{u / v)^Wa{q^)Wi / a{q^) = 5u,vXa , (A.13) 

zeSp ^ k=-i 

since Xa = Wa{z)wi/a{z) is independent of z. By inserting (|A.13I) into the left hand side of 
(IA.12I) we may therefore calculate 

^w^{u/y)wp{y/v)w^{y/w)w^{x/y) = 

= ^^^^Yl Yl '^o.(y/^)^p(y/^)^P/aiy/z)^s/yiz/y')ws{y'/x)w^{y'/w) 
j/GSp zeSp j/'eSp 

= Xa^Yl ^aiv/z)w^{z/u)w^/^{u/v) Ws/^{w/x)Ws{z/w)w^{x/z) . 

The sums over y and y' have been carried out with the help of the star- triangle relation (lA.l II) . 
It remains to recall that Xa^'^s/S'^ / ^) = (^/3/q(^/^))~^ to complete the proof of (IA.12I) . 
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B. Properties of the Q-operator 

B.l Proof of the Baxter equation 

The strategy is similar to [|Ba72llBS90l . Consider 

T(A)-(z|Y(A)|z') = (z|T(A)Y(A)|z'). (B.l) 

The operator T(A) is the difference operator obtained by replacing L^^(\) — )■ L^°(A) in (13.11 ). 
with L^'^(A) obtained from (IB. 21) by replacing u„ and v„ by the corresponding multiplication 
and shift operators u„ and v„ defined in (|2.10l) . 



SG I ~~n \ ~ n ■ n - n ■ n i ■ -n ■ n ■ -n ■ n \ /-r> o\ 



-u ('Q-^^ v"^) A V -A-V"^ 

n V n n " n n I 'n n 'n n 



I 



n 'n n n V^n n n I 



In writing (IB. 21) we have introduced the short-hand notation = ig^K^^and A^ = \/in- Note 
that T(A) acts on the argument z = {zi, . . . , zn) of 1a(z, z'), while it does not act on z'. In 
order to simplify the expression for T(A) we may therefore use a gauge-transformation of the 
form 

Lf(A) = <7„+iL^^(A)(7-\ gn=(^^, °j . (B.3) 
The key point to observe is that 

f Lf (A)2i ■ n(z, zO = X-^ (4 + A„^„uJ (4^, + A„^;iu-i)v-i ■ F,(z, z') 

- A;i (1 + xjj^vr^') (1 + A„t9;iz;+iujv„ ■ r,(z, z') (^-4) 

= 0, 

the last step being an easy consequence of the recursion relations (IA.3I) . (IA.7I) satisfied by the 
functions wa(-2) and w\{z) which appear in the kernel Fa(z, z'). 

Equation (IB. 41 ) implies that 

/ N N \ 

T(A) ■ n(z, z') = I n Lf (A)n + H (A)22 ) " ^(z, z') . (B.5) 

We have 



n=l n=l 



Lf(A)n-r,(z,zO = --^4[C'(^n/< + AAK-A-^(l + ^9„A„^Jz>-i] ■n(z,z') 



' i -t- Z^Z^^iA^I 
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By using the recursion relations (lA.lOl) one may rewrite this as 

Lf (A),i ■ Y,{z, z') = (1/A„, - 1/^„)(1 + q-'XM n ^A^'^- 

where vj^"''^ = w^^^^^^^{zj z'^) «^eA«„/c„(^n4+i)- We may similarly calculate 

= -—{Zn) (An - 1/AA) V„ ■r;,(z,z) 

= - (i/A„ + g/t9J(i - A„^J yj? n ^a'^^- 

It follows that 

N N 
n Lf (A)n ■ >^a(z, z') = n ^(lAn - l/^n)(l + q-'\n^n) " ^^aI^, ^0 , 



n=l ?i=l 
N N 



nLf(A)22-n(z,z') = n^«^n(l/A„ + g/^J(l-AA)-F,,(z,z'). 

n=l ?i=l 

This concludes the proof. 



B.2 Proof of the commutativity 

The key observation to be made is the fact that the operators Y(A) satisfy the exchange relation 

Y(A) ■ (Y(/.*))t = Y(;.).(Y(A*))t. (B.6) 

This is an easy consequence of the exchange relation (IA.12I) . as observed in [IBS90II . Since we 
have A„/ yu„ = A^/ /Um for all n, m = 1, . . . , N we may calculate 

(z|Y(A)(Y(/i*))t|z') = 

N 

ye§N n=l 

N 

= 00 ^ W U)eX„/Kn i^n/ Vn)'^ e^i^/ (l/n/^n)'^eAt„_iK„_i (l/n^n-l)'"^eA„_iK„_i i^n-lVn) 

yg§N n=l 

N 

= '^0 S X\^e^,^/^S^Jyn)w^x^/^Syn/z'n)w,x^_^^^_^{y^z'^_l)W^^^ 

-IKn-l (^n- n/ 

yg§N n=l 

N 

= 5Z n^^Mn/^n(^n/l/n)w^.Mn«n(^n2/n+l)w^.A^K>n/<)^.A*K„(2/n+l<) 
yeSN n=l 

= (z|Y(/.)(Y(A*))t|z'), 
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where 0o = (— as it follows by formula (IA.9I) . The mutual commutativity of 
the operators Q is an easy consequence. Let us furthermore note that (a(A))* = d(A) and 
(T(A))^ = T(A) for A G M. Using (lB?6l) we may calculate 

T(A) ■ Y(A) ■ (Y(/.))t = [a(A)Y(g-iA) + d(A)Y(gA)] ■ (Y(/.))t 

= a(A) Y(/i) ■ (Y(gA))t + d(A)Y(/.) • (Y(g-U))t 

= Y(/x)- [d(A)Y(gA) + a(A)Y(g-U)]^ 

= Y(/.).[T(A).Y(A)]^ 

= Y(/x)-(Y(A))t.T(A), 

which obviously implies [T(A),Q(A)] = 0. 
B.3 Proof of integrability 

In order to prove (13.211) first note that (|3.1 II) allows us to write 

( I Wx{hn) I O = ^k^n\ ^2n I 4 ) Wxiq"^) . (B.7) 
r=-l 

Noting that ( q"^^" \ I Q'^''" ) = ( g^''""^'" | g^''" ) = Sr,k'„-k„ we find that 

{Z^\wx{f2n)\z'r,) = Wx{z^/z'J . (B.8) 

Thanks to this identity and (IA.4I ) it is easy to see that 

N N 

(z I Y(l//te) I z') =YlS,>^,,„w^~2{z^/z'J = (z| JJw«-2(f2„) I z') , 

n=l n=l 

which implies 

N 

Q+(l//s:e) = l[w,-2{f2n)-yl- (B.9) 

n=l 

Similarly note that 

N N 
(z| Y(/t/e) |z') = JJu7«2(2;„4+i) = ( Z | JJ W^2(f2n+l) | z' ) , 
n=l n=l 

which implies 

N 

Q-(/€/e) = Yo ■ nK^(Wi))~' . (B.IO) 



n=l 
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It remains to notice that Yj^ • Yq = Uq to conclude the proof. Indeed, using the notation 

z = {q"^^^, . . . , q'^^''^) and z" = (g^''", . . . , q'^'^'^), we may calculate 

1 ^ 

( z I YJ„ ■ Yo I z") = ^ Yl g"2fc;(fc„+fc„+i)^-2fc;(fcj;+fc;^+i) 

^ {k[,...,k'^)&^ n=l 
N N 

n=l n=l 

= ( z I Uo I z") , 
keeping in mind that we consider the case of odd N. 



C. Asymptotics of Yang-Baxter generators 

From the known form of the Lax operator we derive the following asymptotics for A — +00 
and of the generators of the Yang-Baxter algebras. 

N odd: The leading operators are Bn(A) and Cn(A) with asymptotics: 



(N \ / N (-1)1+° N \ 

n T ^^"^ n ^'^^i"'^" + sub-leading terms, (C.l) 

a=l / \ a=l a=la / 

/ N \ / N (-1)° N \ 

Cn(A) = n - n ^ - n ^-^i"'^"" + ^^b-leadlng terms. (C.2) 

\a=l V V a=l a=l J 

N even: The leading operators are An (A) and Dn(A) with asymptotics: 

(N \ / N (-1)1+" N \ 

n - n + ^'"^ n ^-^^i"'^" + ^^b-leading terms, (C.3) 

a=l V V a=l a=l J 

/ N \ / N (-1)" N \ 

Dn(A) = n - n ^ + A'"* n ^^^^i''^"" + ^^b-leading terms. (C.4) 

\a=l V V a=l ^'^ a=l / 

Note that these asymptotics imply for the SOV-representation of the Yang-Baxter generators the 
following formulae[i] : 

Nodd: 

/ N \ N 

(w^°^)-Mnvi-^)^"Mws°^=n-- (C.5) 



va=l 



"Note that the transformation W^°^ is meant to act as a similarity transformation in the space of the represen- 
tation, i.e. W^'^^ = w^°^/ where w^°^ is a non-trivial operator on space of the states. 
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N even: 

N / N-1 



n^4w^°'')~'0"'w^°'' = hAn^« (c.6) 

a=l V a=l / 

N / N-1 \ 

nea(w^°^)"'eW^°V = [vDll^a]T^, (C.7) 

a=l \ a=l J 

Note that taking the average value of the last two formulae we get for N odd: 

N N 
a=l a=l 

while for N even: 

N-1 N 

ZA={Q)-'l[Z-'l[Xa, ZD = ZA{ey, (C.9) 
a=l a=l 

where (0) is the average value of the charge 6. 



D. Comparison with the Fateev-Zamolodchikov model 

In this appendix we present an explicit comparison between the SG model, studied in this paper, 
and the Fateev-Zamolodchikov lattice model with Zp symmetry nFZ82ll . The Lax operator 
which describes the FZ model has the following expression in terms of the Lax operator of the 
SG model: 

^r(A) = (A)ai. (D.l) 

In the case N(= 2M) even we can construct a map which transforms the transfer matrix of the 
SG model into the one of the FZ model. Let us introduce the unitary operators: 

^n^^n^n = ^nVA = (D.2) 

which in the momentum space play the role of parity operators. Then the unitary operator: 

M 

VTFZ = n (D.3) 

n=l 

has the following action on the Lax operators: 

^FzLi^-aWnpz = icri)'-''Ll^.aii-lY'-^^X) a = 0,1. (D.4) 
so that we get: 



M (A) = a,...M (A)...a, <; ^^^^^^ ^ .,,Qsa^XW,, "^''^ 
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after the flipping ^2n-a — ^ (— 1)^"" ^2n-a of the inhomogeneities. 

In the case N odd the situation is different; the transfer matrices in the two model have different 
spectrum. We use the next two subsections to present an explicit comparison of their spectrum 
in the special case of = 1 and N = 1 . 



D.l Q-spectrum in Sine-Gordon model for = 1 and N = 1 

In this case in the z-representation the operator 2sg(A) is a 3 x 3 matrix^: 

and A± = eA/t^. Now, we observe that 

W,{1) ^ 1, W,iq') = W,iq') = ^ W, 



Ii,ie{l,2,3} 



X + q 

A-1 _ 



Xq — q ^ 



Wa, 



so that in the z-representation: 

esG(A) = 



/ 1 Wa+Wa_ Wa+Wa_ \ 



Wa+Wa_ Wa^ 

V Wa+Wa_ Wa. 



Wa_ 



Then the eigenvalues of Qsg{X) read: 



G)(A) = (Wa^ - Wa J, qr\X) = 1 (1 + Wa^ + Wa_ ± Aa) 

1/2 



(D.6) 

(D.7) 
(D.8) 



(D.9) 



(D.IO) 



with A A = ((Wa+ - 1)2 + 2(Wa+ - 1)Wa_ + (1 + Swl^)wl_^ and clearly (2sg(A) has 
simple spectrum for all the values of the local parameter k, E C 



D.2 Q-spectrum in Fateev-Zamolodchikov model for q^ = 1 and N = 1 

In this case in the z-representation the operator 2fz(A) is a 3 x 3 matri5<0: 

eFz(A) = ||(^ = q'^'-'^\Q,z{X)W = g2(^-^)) = l^A(g'(^-^'^)WA(g2(^-^-))||,,e{i,2.3}, (D.ll) 



'^Note that to make more simple the comparison with the Q-operator of the FZ model, here we have considered 
for QscW the operator Y(A) defined in ( 13.101 ) just with a different normalization. 
'^Here we have rewritten in our notation the (5.12) of I1BS901 for fc = 0. 
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(D.12) 



(D.13) 
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